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1 Introduction 



Two-dimensional tethered surfaces, which model polymerized flexible membranes, offer 
interesting problems of statistical mechanics (for a general introduction see [f]). One 
problem, which is now reasonnably well understood, is the influence of bending rigidity: 
for high bending rigidity (or low temperature), such a membrane is flat, with fractal 
dimension dj = 2, while for low bending rigidity (or high temperature) it is crumpled, 
with infinite df. The existence of a crumpling transition, separating these two phases, 
has been established by numerical simulations and by renormalization group calculations 
for "phantom surfaces", where self-avoidance is ignored. 

The effect of self-avoidance and its interplay with bending rigidity is more difficult 
to study and it is not so well understood. Numerical simulations of self-avoiding flexible 
polymerized surfaces in 3-dimensional space favour the idea that self-avoidance has drastic 
consequences and flattens the surfaces at any temperature [2, 3, 4]. There is however no 
fully convincing analytical argument for such a behavior. Another question is how the 
behavior of self-avoiding surfaces depends on the dimension of bulk space or on the details 
of the contact interaction. 

The standard model for theoretical studies of self-avoiding surfaces has been first 
discussed in [5, 6] and is inspired by the Edwards model [7] for polymers (for a general 
presentation see e.g. [8]): it consists in an extension of this model from a line (the 
polymer) to D-dimensional manifolds. The case D = 2 corresponds to surfaces. The 
internal points of the manifold, which belong to the nodes of a D-dimensional network, 
are labelled by continuous coordinates x £ IR . The position of these points in the 
external o?- dimensional bulk space is described by the vector field x — > r(x ) £ U d . The 
continuum Hamiltonian is 



The first term is the Gaussian elastic term, which describes the crumpled phase of "phan- 
tom" surfaces. The second term is a "weak" self-avoidance 2-body <5- potential, which 
models the contact interaction in bulk space, t > is the coupling constant. 

Dimensional analysis shows that the contact interaction is relevant at large distances 
if d < o?*, where d* = AD/ (2 — D) is the upper critical dimension. As in the case of 
polymers, it is natural to perform an e-expansion about d* to evaluate scaling properties 
such as the fractal dimension df. Since for D = 2, the upper critical dimension d* = oo, 
it is in fact better to perform both an expansion in d and D } starting from (o?*, D*) with 
D* < 2, aiming e.g. at d = 3, D = 2. 

The first calculations for the model (1.1) [5, 6] used the direct renormalization method, 
adopted from polymer theory [9, 10]. It has been developed by several authors to perform 
calculations at 1-loop order [11, 12, 13]. In this method the theory is reexpressed in terms 
of dimensionless physical quantities, determined for finite surfaces with internal extent L. 
This length L provides a renormalization scale and allows to calculate the renormalization 
group flow for the model at 1-loop order. The method has been checked to be valid at 
1-loop order [14]. 

Recently a different and more general formalism has been introduced by B. Duplantier, 
E. Guitter and one of the authors in [15]. (It partly relies on previous studies of mem- 
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branes interacting with a fixed element [17, 18, 19]). Although the interaction term in the 
Hamiltonian (1.1) is a non-local and singular function of the field r(x), it is shown that 
the short distance behavior of the model can be encoded in a multilocal operator product 
expansion (hereafter abbreviated as MOPE), which generalizes the Wilson operator prod- 
uct expansion valid for local field theories. This allows a systematic analysis of the short 
distance ultra-violet (UV) singularities of the model and shows that the theory (1.1) is 
renormalizable at the critical dimension d*. This means that the model is rendered UV 
finite in perturbation theory at d* by a renormalization of the coupling t and the field r 
("wave-function" renormalization). In parallel with the derivation of the renormalization 
group equations for the $ 4 -theory in dimension d = 4 — e < 4, which gives the scaling laws 
for a large class of critical phenomena, renormalization group equations for the model of 
self-avoiding tethered surfaces are derived in a similar expansion for d < d*. The 1-loop 
results obtained through this method confirm the previous calculations of [5, 6]. The 
consistency of the direct renormalization method follows from the validity of finite size 
scaling laws for finite manifolds, also established in [15]. 

In this paper we apply the renormalization group approach of [15] to a different prob- 
lem: that of tethered surfaces at the tricritical point, the so-called ©-point. Our motiva- 
tion is twofold: 

• Physical: the study of the ©-point for membranes is physically interesting in its own: 
For polymers it exists due to a competition between 2-body attractive interactions 
(like long-range Van der Waals forces) and hard-core repulsive interactions. At high 
temperature the repulsive interactions dominate and the polymer is swollen. At 
low temperature attractive interactions dominate and the polymer is in a collapsed 
compact state. For a single long polymer, the transition between these two states 
occurs at the ©-point. This point represents a different multi-critical state for the 
polymer [20] (we refer to [8] for a general presentation). One expects a similar 
transition to occur for membranes and an interesting fact was first pointed out in 
[15]: for polymers and for membranes with internal dimension D small enough, one 
expects an effective 3-body repulsive interaction to be relevant to describe the 0- 
point close to the upper critical dimension d c = 3D/{2 — D). For higher D } it is a 
modified 2-body interaction, repulsive at short range, but attractive at larger range, 
which is relevant to describe the ©-point close to the upper-critical dimension, now 
given by d' c = 2(3D — 2)/(2 — D). The crossover between the two interactions occurs 
at D = 4/3, d = 6. While the modified 2-body interaction leads, at 1-loop order, to 
calculations for the critical exponents which are analytically computable and quite 
similar to those for self-avoiding membranes, the 3-body interaction has not been 
considered up to now - except of course for polymers (D = 1) which have been 
extensively studied - and it is not known which theory should describe "physical 
membranes" with D = 2 and d = 3. As we shall show here, the interplay between 
the 2-body and the 3-body interaction can be studied via a "double e-expansion" 
around the critical point D = 4/3, d = 6. 

• Mathematical: 1-loop calculations for the model of membranes with the repulsive 3- 
body interaction are already non-trivial and the first order term of the e-expansion 
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appears not to be computable analytically, i.e. cannot be expressed in terms of 
standard special functions except for D = 1. In fact, at 1-loop order, one encounters 
problems similar to those occuring in the evaluation of 2-loop corrections for self- 
avoiding surfaces but without the additional difficulty of double poles. Thus, this 
model can be considered as a (not so enjoyable) toy-model to develop analytical 
as well as numerical technics which should apply to higher order calculations for 
self-avoiding surfaces. 

This paper is organized as follows: In section 2 we define the model with the 3-body 
interaction, recall how its perturbative expansion is obtained and how the short distance 
UV divergences are organized according to the MOPE. We then show explicitly how these 
divergences can be subtracted in order to construct a renormalized theory and how the 
scaling laws are obtained. 

Sections 3-6 are devoted to the explicit calculation of the 1-loop counterterms for the 
3-body interaction, valid a priori for D < 4/3. 

The counterterm associated to the "wave-function" renormalization, i.e. to the renor- 
malization of the elastic energy term in the Hamiltonian, is treated in full details in 
section 3. The integral representation of the counterterm in terms of the MOPE coef- 
ficient is derived. Various technical problems are discussed: extraction of the singular 
part (residue of the pole in e), subtraction of divergences associated to the fine-tuning of 
the 2-body interaction needed to reach the ©-point and the definition of the measure in 
non-integer dimension D. Using these methods, the counterterm is evaluated numerically 
in the range 1 < D < 4/3. Finally the result for D — > 1 is compared with the value for 
D = 1, already known analytically [21]. 

In section 4 we introduce a useful change of variables for the integral representation 
of the counterterms, based on conformal transformations. 

In section 5 the first two diagrams for the coupling constant counterterm are discussed. 
The first one is calculated analytically. The second is evaluated by methods similar to 
those used in section 3 for the wave-function counterterm. 

Section 6 is devoted to the evaluation of the last diagram needed for the coupling 
constant counterterm. It can only be evaluated numerically. A full use of the technical 
tricks developped previously (analytic continuation of the integration measure, conformal 
mappings) is required as well as the implementation of an original adaptative Monte-Carlo 
integration method. 

After this rather technical part we gather in section 7 the various counterterms for the 
model with 3-body interaction. We further give the results for the anomalous correction 
at 1-loop order to the exponent v in the range 1 < D < 4/3. This exponent is related to 
the fractal dimension of the membrane at the ©-point. 

In section 8 we study the crossover between the 3-body interaction and the modified 
2-body interaction. First we show that the 1-loop corrections obtained previously from 
the 3-body interaction have a smooth limit for D — > 4/3. Then we recall the 1-loop 
corrections obtained from the modified 2-body interaction (valid for D > 4/3) and show 
that the limit D — > 4/3 exists and is close to but different from the previous one. Finally 
we show that the interplay between the two interactions can be studied for D and d 
close to the critical values D = 4/3, d = 6, by a new "double e-expansion" . At 1-loop 
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Figure 1.1: Critical dimensions for various operators (solid lines) and phase diagram. The phase separa- 
trices are the fat lines. 



order this expansion is analytically computable and the renormalization group flow can 
be studied explicitly. As a result we show that depending on D and d, the ©-point is 
described either by: (a) a Gaussian fixed point, (b) the 3-body repulsive interaction, (c) 
the modified 2-body interaction. The three corresponding domains in the 2-dimensional 
(d,D) plane are depicted in figure 1.1. The fat lines are the separatices between these 
domains. The fat dashed line is a linear extrapolation of the 1-loop result. This line 
separates the domains (b) and (c). As discussed in section 8, this result indicates that 
the modified 2-body interaction should be relevant to describe 2-dimensional membranes 
at the ©-point, independently of the dimension d. 

The results are summarized in the conclusions. More technical points are discussed in 
the appendices. 

Appendix A treats problems associated with the finite part prescription which we use 
to subtract the relevant UV divergences. It is shown that different prescriptions may be 
adopted but lead to the same 1-loop results. 

In appendix B diagrams not calculated in the main text are given. 

Appendix C briefly discusses the anomalous dimension of the 2-body self-avoiding 
interaction. This allows to describe the model in the neighborhood of the tricritical 
point. 
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2 The Model 



2.1 The 3-body Hamiltonian 

The field r(x), r £ IR d , x £ 1R D describes the conhguration of the D- dimensional poly- 
merized flexible membrane in o?- dimensional space. The Hamiltonian for this membrane 
with 3-body repulsive interaction is 

^°H= / Uvr(x)) 2 +g III S d (r(x)-r(y))S d (r(x)-r(z)) (2.1) 

and J x = J d D x. The calculations of physical observables are performed as an expansion in 
g and an analytical continuation in the internal dimension D and the external dimension 
d (dimensional regularization), along the line of [15]. Dimensional regularization allows to 
deal with the short distance (ultraviolet) divergences, which appear as poles in the complex 
D or d planes. Otherwise a physical short distance regulator has to be used. One can for 
instance add a term proportional to the curvature to the Hamiltonian, which in the spirit 
of the Pauli- Villars regularization contributes higher derivatives to the elastic energy term 
(Vr) 2 , thus modifying the free propagator at short distances. Large distance (infrared) 
divergences also occur. They can be cured by using an IR regulator, for instance by 
considering a finite membrane. However the calculations become technically more difficult 
and one has to keep track of curvature, boundary and finite-size effects. Alternatively 
observables invariant under global translations (r(x) — > r(x) + r ) in bulk space may 
be considered, since these observables are expected to be IR finite even for an infinite 
membrane. For more details cf. [16]. 

Dimensional analysis shows that the dimension of r and of the coupling constant g are 
(in internal momentum units such that [x] = —1) 

[r] = -u = , [g]= £ = W- 2vd (2.2) 

The interaction is relevant in the sense of Wilson if e > and the perturbation theory is 
expected to be UV finite up to subtractions associated with relevant perturbation terms, 
which we shall discuss later. The interaction is irrelevant if e < 0. The short distance 
divergences will occur as poles in e at e = 0. As for standard Landau- Ginzburg- Wilson 
$ 4 models, which describe critical phenomena in the e = 4 — D expansion, these poles 
have to be subtracted in order to define a renormalized held theory UV finite at e = 0. 
This theory will give the scaling behavior of the model for e > 0. 

For clarity we shall graphically represent the different interaction terms which have to 
be considered. The local operators are 



1 = • (2.3) 
Vr(;c)) 2 = + . (2.4) 



1 - ^2 
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The bi-local operators are 



6 d {r(x) - r(y j) = • • (2.5) 

(-A r )6 d (r(x) - r(y)) = • — h — •. (2.6) 
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The tri-local operator in the Hamiltonian (2.1) is 



8 d (r(x) — r(y))8 d (r(x) — r(yj) 



{2.7) 



The hrst terms of the perturbative expansion in g of the expectation value of an 
observable O are 



(O)o 

1 2 

+ - Q 



(O 
(O 




)o + 



f2.8) 



where (O)o will denote from now on the expectation value of an operator O for the free 
theory with g = 

(O) = J [ U . (2.9) 

J V[r]e-Jx + 

The perturbation expansion (2.8) will suffer from short distance singularities, which occur 
when subsets of points coalesce. This gives rise to the renormalization discussed in the 
following. Examples of IR finite observables are provided by "neutral" products of vertex 
operators 

N 



N 



o 



U 

a=i 



. 



(2.10) 



a = l 



To compute the (. . .) in equation (2.8), one writes the ^-functions of the interaction 
operator as the Fourier transform of vertex operators, performs the free average, then 
inverses the Fourier transformation and ends up with an integral over the positions of the 
internal points belonging to the interaction operators. The integrand is a function of all 
the distances between the internal points and the external points (the x a } s). We will show 
in 3.1 how this works. 



2.2 The MOPE: mixing of 3-, 2- and 1-body operators 

Short distance singularities may arise when the distances between internal points vanish. 
In fact, for observables of the form (2.10), no additional singularity occurs when distances 
between internal and external points vanish as long as the distances between external 
points stay non-zero. As shown in [15, 16], the short distance behavior of expectation 
values of operators in the free theory is given by a multilocal operator product expan- 
sion. This implies that the short distance divergences can be absorbed by adding to the 
Hamiltonian (2.1) counterterms proportional to multi-local operators. 

At hrst order in g, these operators are generated by contracting points in the single 
3-body operator. If the three points (x\ } x 2} x 3 ) are contracted towards their center of 
mass o, 1-body operators are obtained: 




= A({x 1 ,x 2 ,x 3 }) • + B({x 1 ,x 2 ,x 3 }) + + ... (2-11) 
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where {x\ } x 2} x 3 } means the relative distances between the 3 points X{. By power counting 
the coefficients A and B are homogeneous functions of the difference between the positions 
(xi, x 2} x 3 ) of the points, with degree given respectively by 

deg(A) = -2vd , deg(B) = D — 2vd . (2.12) 

In equation (2.11) we made use of the equation of motion of the theory and neglected 
total derivatives such as V(rVr), since they do not give UV divergences. The dots . . . 
represent operators of higher dimension. 

If only two points are contracted (e.g. Xi 7 x 2 ) } 2-body operators are generated: 




- = C({x u x 2 }) • . + D({x u x 2 }) • h . + ... (2.13) 



with 



deg(C) 



'd 



deg(D) 



d)i 



(2.14) 



At order g 2 } one can contract subsets of points of a pair of 3-body operators. Since 
we are interested in evaluating anomalous dimensions at 1-loop order, we only need to 
know the divergences which arise when two 3-body operators coalesce into a single 3-body 
operator. Three contractions are possible: 



£({^i,2/i}, {^2,2/2}, {^3,2/3}) 



+ 



(2.15) 




F({xi, 2/i, 2/2}, {^2,2/s}) 



£({£2,^3,2/2, 2/3}) 




with 



+ 



+ 



deg(E) = deg(F) = deg(G) 



-2vd 



(2.16) 



(2.17) 



(2.18) 



The coefficients A in equation (2.11) and C in equation (2.13) give strong short distance 
divergences. We call a divergence strong or relevant if it is not integrable and thus has to 
be treated by a finite part prescription. The relevant divergence here is expected, since 



the corresponding operators 1 and 6 a 



r x ; 



r(y)) are relevant for e = 0. The divergence 



proportional to the unity operator 1 is just a "vacuum energy" term. It does not occur for 
expectation values of physical observables, but will be present in the partition function 

of finite membranes. The divergence proportional to the 2-body operator • • has 

to be cancelled by adding a 2-body counterterm 



AH°\ 



At 



8 d (r(x) — r(y)) 



(2.19) 



x Jy 
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and by fine tuning At so that the renormalized 2-body interaction vanishes. This sit- 
uation is known from the ©-point of polymers: the fine tuning of the 2-body interac- 
tion is required in order to reach the ©-point which separates the swollen phase, where 
self- avoidance is relevant, from the collapsed phase, where attraction and short-distance 
repulsion dominate. It also arises in standard scalar held theories: Quadratic divergences, 
associated to a mass renormalization, have to be subtracted in order to stay at the critical 
point. 

When using dimensional regularization to define the theory, these operators do not 
give rise to logarithmic divergences at e = and thus can be subtracted unambiguously 
by a finite part prescription. This amounts to analytically continue beyond the poles 
caused by these operators. With this prescription, the renormalized coupling of the 2- 
body interaction is automaticaly zero if the bare coupling t is set to zero, i.e. if one starts 
with the 3-body Hamiltonian (2.1). 

The subleading term D({x\ } x 2 }) in equation (2.13), proportional to the modified 2- 
body interaction (2.6), gives a divergence for e = if the internal dimension D of the 
membrane is larger than or equal to 4/3. This reflects the fact, first outlined in [15], 
that at the ©-point, the modified 2-body interaction is more relevant than the 3-body 
interaction if D > 4/3, while for D < 4/3, including the case of polymers, the 3-body 
interaction is the most relevant one. The Hamiltonian (2.1) thus describes the ©-point 
for D < 4/3 and e small. For D > 4/3 it describes a multicritical point reached by fine 
tuning both the 2-body and the modified 2-body couplings. A more serious investigation 
of the relative relevance of these two operators for finite e requires a study of the model 
with both couplings, 

Kb = i + + 9 I I I , + b j j , (2.20) 

J x J x Jy J z \^ J x Jy 

around the point D = 4/3 and d = 6. This topic is discussed in section 8. 

2.3 UV divergences and 1-loop renormalization 

Let us now concentrate on the Hamiltonian (2.1) and on the divergences at e = 0. The 
MOPE structure of the UV singularities implies that the theory can be made UV finite 
for e — > by considering the Hamiltonian 

K[r] = Z j Uvrix)) 2 + gZ g ^ j j j 6 d (r(x) - r(y))6 d (r(x) - r(z)) , (2.21) 

where // is a renormalization momentum scale, r and g are the renormalized field and 
the renormalized coupling constant. As in [15, 16], the counterterms Z (wave-function 
renormalization) and Z g (coupling constant renormalization) subtract the poles at e = 0. 

The counterterms are evaluated as follows. The divergence proportional to the oper- 
ator -j- occurs because of the integration over a global length scale in the MOPE (2.11). 
For instance, if we integrate over the three points (xi, x 2} x 3 ) in a domain T>l(o) of size L 
around their center of mass o 

^l(o) = {^i + x 2 + x 3 = 3o and \x a — Xb\ < L } a, b = 1, 2, 3} , (2.22) 
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using equation (2.11) gives a pole of the form 



L £ 

e 



+ } + + 0(e°) 



(2.23) 



with the residue determined by the coefficient B of the MOPE and abbreviated graphically 

as 



+ ) = //* 1+ * 2+ * 3 =o B({x u x 2 ,x 3 }) 

sup(\x a -x b \) = L 



(2.24) 



e=0 



We will explain that in more detail in section 3.3. 

The other residues come from a similar integration at a "typical" distance L between 
the points in the clusters of the coefficients E, F and G in equations (2.15), (2.16) and 
(2.17). They are abbreviated similarly as 



up of distances 
n clusters = L 



£({^i,2/i}, {^2,2/2}, {^3,2/3}) 



e=0 



(2.25) 



up of distances 
n clusters = L 




^({^1,2/1,2/2}, {^2,2/3}) 



e=0 



(2.26) 



(2.27) 



e=0 



/ ... t G({x 2 ,x 3 , j/ 2 , y 3 }) 

'sup of distances v ^ J 7 

in cluster = L 

With these notations and using equation (2.8), the counterterms which make the 
theory hnite at 1-loop order are 



Z = 1 + g- +0(g 2 

e 



+ 



(2.28) 



1 + g- +0(g 2 

e 




+ 18 



9 

+ 2 




(2.29) 



The coefficients of the counterterms, a and 6, depend of course on the point (D C} d c ) on 
the curve e = 3D — 2vd = 0, where the renormalization is performed. 

The picture associated with this renormalization prescription at 1-loop order is to cut 
out domains of size L = ^ around collapsing points in the perturbation expansion (2.8). 
The freedom in the choice of a renormalization scale // is represented as the freedom to 
choose the size L = - of these domains. 
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2.4 RG equations and scaling relations 



From the existence of counterterms and of a UV finite perturbation theory for e > one 
deduces renormalization group (Callan-Symanzik) equations for the renormalized theory 
and scaling laws for the model for e > in the standard way. The renormalized Hamilto- 
nian (2.21) can be rewritten Hamiltonian (2.1) through the change of variables 

to bare held r (x) and bare coupling constant g 



r (x) 



Z 1/2 r(x) 



go 



gZ Z g p £ 



(2.30) 



The Callan-Symanzik equations, which give the scale dependence of the renormalized 
theory, are obtained via the g dependence of the renormalized couplings keeping the bare 
couplings fixed. One thus obtains the renormalization group /3-function for the coupling 
constant 

d 

g = -eg + (ad+b)g 2 + 0(g 3 ' 



P B (9) = ^ 



(2.31) 



90 



and the scaling dimension v of the held r 



1 d_ 



InZ 



90 



-7f~ + \ a 9 + °{9 2 ) 



(2.32) 



As we shall see, both a and b are positive and the renormalization group flow has an IR 
stable fixed point for positive g = g*. = v(g*) 
the membrane at the 0-point via 

D 



is related to the fractal dimension dp of 



d% = — 



(2.33) 



2.5 A change in normalizations 



Before describing the details of the calculations for the counterterms, let us change the 
normalizations used in the definition of the theory. This is done for technical purpose 
only, but appears convenient to avoid a lot of factors tt and T functions in intermediate 
expressions. We rewrite the bare Hamiltonian as 



n°Jr] 



2-D)- 1 / + + g 




with as before 



1 / \2 



but with a modified integration measure 

1 



s 



D 



D 



T(D/2) 



(2.34) 



(2.35) 



(2.36) 
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being the volume of the (D — l)-dimensional unit sphere. With this normalization the 
free propagator becomes simply 



(r(x)r(y)) = - \x - y\ 2v . 
Similarly, we normalize the measure in Fourier space ]R d as 

J = K- d ' 2 J d d p 

and the (^-distribution in ]R d is modified to 

6 d (r-r') = (ivf 2 6 d (r -r r ) = 

so that the 2-and 3-body operators, which have to be considered, are 

• • = S d (r(x) — r(y)) , 



• = (-A r )6 d (r(x)-r(y)) , 

= 8 d (r(x) — r(y))8 d (r(x) — r(zj) . 

The only change in the expressions for the 1-loop counterterms (2.28) and (2.29) i 



(2.37) 
(2.38) 

(2.39) 

(2.40) 
(2.41) 

(2.42) 



(2-D) 



+ 



(2.43) 



3 Elastic Term Renormalization 

3.1 Explicit evaluation of the 3-point MOPE coefficient 

Let us hrst derive explicitly the MOPE (2.11) for the 3-body operator (2.42). The 3-body 
operator has as Fourier integral 



ipr(x) .. e iqr(y) .. p —i(p + q)r(z) . 



:: e 



(3.1) 



v J <i 



: : denotes the usual normal product for free fields. We use the OPE for products of vertex 
operators 

. e ipr(x) . . e iqr(y) . . e ~i(p + q)r(z) . 

= . e ipr(x) e iqr(y) e -i(p + q)r(z) . ^ipr(x) + iqr(y) - i(p + 9)^(^)^(3,2) 

and since r is a free held we have explicitly, using (2.37) 
l^ipr(x) + iqr(y) - i(p + q)r(z)^ 

= e -\{(p( r i x ) - r i z )) + i{ r {y) - r i z ))f)o 



p 2 \x — z\ 2v + (pq) (\x — z\ 2v + \y — z\ 2v — \x — y\ 2v ) — q 2 \y — z 



\2v 



(3.3) 
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a 



y 



Figure 3.1: The triangle A(a,b,c) 



Let o be the center-of-mass of the 3 points x, y and z and let us denote 

x = o + x , y = o + y , 2: = o + 5 , 5: + ?/ + 5 = (3-4) 
The MOPE is obtained by performing a Taylor expansion around the center-of-mass o 
. e ipr(x) e iqr(y) e -i(p + q)r(z) . = e xd x Jjdy e zd z . e ipr(x) e iqr(y) & -i(p + ?)r(z) 



Expanding in powers of d } we find up to order 2 in d 
e xd x Jjdy e zd z . e ipr(x) + iqr(y) - i(p + q)r(z) 



x—y—z—o 

(3.5) 



(3.6) 



x=y=z=o 



-\x - z\ 2 p 2 + 2(y - z)(z - x) (pq) - \y - z\ 2 q 2 1 2 
= Ho) + -:(Vr) (o):+... 

The factor 1/dD comes from the contractions between internal space indices (in the partial 
derivatives d x and x) and between external space indices (in r and the external momenta 
p). 

Using (3.3) and (3.6) and performing the explicit integration over p and q } the term of 
order will give the coefficient A in (2.11), whereas the term of order 2 yields B. Denoting 
by a, 6, c the respective distances between the three points x, y and z (see figure 3.1) 

a = \y — x\ , b = \z — y\ , c = |x — z| , (3-7) 

we find 

A({x,y,z}) = I f^Pr(x) + iqr(y)-i(p + q)r(z) )o 

Jp J q 



and 



B({x,y,z}) 



(a" + b v + c i/ )(a i/ + ^ - c i/ )(6 i/ + c v - a v )(c v + a" - b u ) 
— \x — z\ 2 p 2 + 2(y — z)(z — x) (pq) — \y — z\ 2 q 2 



(3.8) 



v J <i 



dD 



l^ipr(x) + iqr(y) - i(p + q)r(z)^ 



(3.9) 



2 d (c 2 + a 2 - b 2 )b 2v + (a 2 + fe 2 - c 2 )c 2 » + (fe 2 + c 2 - a 2 )a 2 ^ 
D \a v + b l/ + c v )(a v + b v - c v )(b v + c v - a v )(c v + a v - b v )\ 1+ '' 
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3.2 MOPE and relevant subdivergences 



If two out of the three points in I / are contracted, a dominant subdivergence occurs. 
For instance, if x — > y, the coefficient B given by equation (3.9) behaves as 

B({x,y,z}) ~ -^\x-y\- vd \x-z\ D -» d (3.10) 

As explained in section 2.2, this follows from the MOPE (2.13). Using similar techniques 
as for calculating the coefficients of the MOPE (2.11), the leading coefficient C in equation 
(2.13) is found to be 

C({x,y}) = \x-y\-" d (3.11) 

To subtract this divergence with a finite part prescription, we must add a 2-body coun- 
terterm (2.19) which is at first order in g 

or equivalently we can redefine the 3-body operator (2.42) as 
= 8 d (r(x) — r(y))8 d (r(x) — r(zj) 



*|-"F(r(*) - r(,)) !(;) + P 8 ™-"^ (3.13) 
2 [ of x,y, z J 

With this subtraction prescription the coefficient B in the MOPE (2.11), given originally 
by equation (3.9), becomes 

B({x,y,z}) = (3.14) 
2 d (c 2 + a 2 - b 2 )b 2v + (a 2 + b 2 - <?)<?" + {b 2 + c 2 - a 2 )a 2v 



D \a v + b v + c v ){a v + b v - c v ){b v + c v - a v ){c v + a v - b^ * ' 

f jyD — vd _|_ ^D — vd ^D — vd _|_ ^D — vd ^D — vd _|_ j^D — vd ^ 
"I - . -i-i i ..j "I - 7.. j "I - 



4£> [ a 1 ''' ' 6^ ' J 

This can be checked by an explicit calculation along the lines of section 3.1. 



3.3 IR regularization and extraction of the residue 

As discussed in [18], integrals over points x, y } . . . on the membrane are defined in non- 
integer dimension D by switching to distance variables, that is by integrating over the 
relative Euclidean distances \x — y\ } ... between these points. Through this change of 
variables the usual measure d D x A dPy A . . . becomes a distribution over the distance 
space, which depends analytically on the dimension D. This distribution will be shown 
more explicitly in section 3.4. 
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We already discussed that UV- divergences have to be treated by dimensional regu- 
larization. The problem arising in this context is that our Hamiltonian has no intrinsic 
length scale. Choosing the dimensions D and d so that e > in order to make Feynman 
diagrams UV-hnite thus necessarily involves IR-divergences. To extract the counterterms 
an IR-regulator has to be introduced. This regulator will set the renormalization scale. 
We want to extract the UV divergence associated to the global contraction of subsets 
of points. An example of such a contraction, involving four subsets indicated by dotted 
lines (one subset is reduced to one point), is depicted in figure 3.2. According to [15] this 




Figure 3.2: A contraction of points into 4 subsets and the corresponding MOPE coefficient for 4-body 
operators. 

divergence is proportional to the marginal multilocal operators appearing in the MOPE 
for this contraction. The divergence is obtained by integrating the corresponding MOPE 
coefficients for small distances inside each subset and by keeping the logarithmically di- 
vergent part. By construction, the MOPE coefficient C depends only on the relative 
distances between points inside each subset and not on the relative distance between the 
different subsets. (This is why in figure 3.2 the coefficient C is noted as an independant 
function of each subset; in fact the last subset, reduced to one point, does not play any 
role in C .) 

As schematically discussed in section 2.3, we choose the following IR-regularization 
prescription: The distances in any contracted subset have to be smaller than L, the 
IR-cutoff. Compared to other prescriptions like considering closed membranes with the 
topology of a hypercube or a hypersphere of diameter L, this prescription has a great 
calculatory advantage. The integrands and the MOPE coefficients are not modified and 
the interesting pole term can be extracted simply by using homogenity. 



Let us apply this regularization prescription to the graph 
The MOPE (2.11) gives with the conventions of figure 3.2: 



that we consider here. 



(3.15) 



When integrating B over the three points x, y } z with fixed bary center o and with the IR 
regulator L, i.e. in the domain 



^l(o) = {x + y + z = 3o and \x — y\ } \y — z\, \z — x\ < L} 



(3.16) 
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we obtain the UV divergent term, which is denoted 



+ 



B[ ... 



L Jv l{o) 

Since the integrand (3.14) is homogenous, we deduce that 



(3.17) 



+ 



c e 



(3.18) 



For e = 3D — 2vd > the integral (3.17) is UV and IR convergent and the residue c(e) is 
analytic in e. We are interested in the residue c(0). Applying Z^- on both sides yields: 



L 



d 



(3.19) 



Acting on the r.h.s. of equation (3.17) the operation L-^ extracts the boundary dT>i,{o) 
of T>l(o). There the largest distance between the three points equals exactly L. 



c(0) 



B( 



dv L (o) 



(3.20) 



e=0 



The normalization introduced in equation (2.36) was chosen in order to eliminate the 
additional factor Sd one might expect here. 

We can consider separately the contribution of each so-called sector, where a = \x — y\, 
b = \y — z\ orc = \z — x\ respectively is the largest distance (the set where two distances 
have exactly the same length is a set of measure zero and thus can be neglected). By 
symmetry each sector gives the same contribution. For c(e) the prescription yields: 



ce 



2 d 
~D 



1 



J L=a>6,c 

(c 2 + a 2 - b 2 )b 2v + (a 2 + b 2 - <?)<?" + {b 2 + c 2 



(3.21) 



a a 



(a v + b v + c v )(a v + b v - c v )(b v + c v - a v )(c v + a v - ¥ 

D-vd _|_ yj-vd 1 



1 + ^ 



jyD — vd _|_ ^D — vd 



,vd 



C D- V d + a D-vd a . 



ivd 



The hnal expression for the residue c(0) is obtained by replacing d by d c = 3-D/ (2 — -D) 
in equation (3.21). That this replacement is justified will be discussed in appendix A. 



3.4 Analytic continuation of the measure 

We now define the explicit form for the integration measure in non-integer dimension _D, 
that will be used in the calculations. For other but equivalent formulations we refer to 
[18]. 
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The general problem is to integrate some function / invariant under translations and 
Euclidean rotations over all configurations of N points x X} . . . , xjv imbedded in D dimen- 
sions. This implies that / depends only on the N(N — l)/2 relative distances \xi — Xj\ 
between these points. In the following the integral over the center of mass is therefore 
always excluded. In order to be able to define the integration, let us take D > N — 1 and 
integer. For i < N denote by = X{ — x^ the z'th distance- vector and by y? its a'th 
component (a = 1, . . . , D). 

The integral over y x is simple: Using rotation invariance, we fix y x to have only the 
a = 1 component non-zero. The measure becomes 

/rco 
d% = S D J o dy{ (ylf- 1 , lh = (yj, 0, . . . , 0) (3.22) 

where Sd is the volume of the unit-sphere in IR D , defined in (2.36). 

We now fix y 2 to have only a = 1 and a = 2 as non-zero components. The integral 
over y 2 consists of the integration along the direction fixed by y x and the integration in 
the orthogonal space IR D_1 : 

/yoo /'CO 
d% = S D -! / dyl / dy 2 2 (y 2 2 ) D - 2 , 2/2 = (yl yl 0, . . . , 0) (3.23) 
J -co JO 

For the j-th point, one proceeds recursively to integrate first over the hyperplane defined 
by j/i, ... , yj-i and then the orthogonal complement: 

//'CO /'CO 
d% = S D - j+1 n / '/.'/^ ( '/.'/ ijJ^f- 3 , // i//' //; •" 0) (3.24) 

The final result for an integral over all configurations of N points is 

/■N-l N-l roo . . \ 

/ n ^ = sdSd-! . . . s D . N+2 n n / ^ m wf~ 3 ^ 

J j = l j = l \a = l J -°° J ° J 

This provides a well defined analytic continuation of the measure to dimensions D non- 
integer, even for D < N — 1. It is equivalent to the measures defined in [18]. 

As we first want to contract 3 points with coordinates (xi, x 2} x 3 ) } let us look what the 
general expression results in: By translation invariance we fix one point, x 3 and integrate 
only over y = x\ — x 3 and discussed above. We thus have to evaluate an 

integral of the form 



d y d z f(y,z) = S D S D -i dy 1 y 1 / dz t / dz 2 z 2 f(y 1 ,z 1 ,z 2 ) (3.26) 

JO J-co JO 

This expression is well defined and integrable for D > 1. For — > 1 it should reduce to 
the integral over a line. Let us take D = 1 + e then we get 

/•CO /'CO /'CO 

2e / dyi / dz x \ dz 2 z e 2 ~ x f(y 1 , z t , z 2 ) , (3.27) 

JO J-00 JO 

where subleading terms in e from the expansion of SdSd-i are neglected. As ez 2 ~ x is a 
representation of the ^-distribution for e — > 0, (3.27) reduces in this limit to 



2 ^ dy x j dz 1 f(y 1 ,z 1 ,0)= j dy x \ d Zl f(y u Zl ,0) (3.28) 



as expected. 
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3.5 Explicit integration 

From (3.21) the explicit expression for the pole term proportional to ~ e L e of the diagram 
contributing to the wave function renormalization is: 



+ 



Id 



( a 2" + <?» _ b 2v )b 2 + (c 2 ^ + b 2v - a 2v )a 2 + {b 2v + a 2v - c 2v )c 2 

l=a>c>b L J 

4 \ <ic/2 + l 

(a i/ + c i/ + 6 i/ )(a i/ + c i/ -6 i/ )(a i/ -c i/ + 6 i/ )(c i/ + 6 i/ -a i/ ) J 
-a- udc {b D - udc + c D - udc ) - b' 
-c- vi *{a D - vi * + b D - vdc ) (3.29) 



-vd c ( „D — vd c i „D — vd. 



[a 



The factor 6 is due to the explicit ordering of the distances a, b and c. This expression 
has to be integrated numerically. From (3.23) and with the change of normalization in 
(2.36) we know that the measure is 



/(a, 6, c) = / db x / db 2 b 2 D ~ 2 f(a,b,c) 

OD J — oo Jo 



a = 1 



6 = A + b 2 , c=^(l-6 1 ) 2 + & 



2 

2 5 



(3.30) 



restricted to the domain, where b < c < a: 




Figure 3.3: The half-sector used for the numerical integration 



For 1 < D < 4/3, the expression (3.29) is integrable everywhere. For D < 1, the measure 
has a non-integrable singularity. For D > 4/3, an additional non-integrable singularity 
appears for small b. 

Various numerical problems exist: The hrst is due to the integrable singularitiy for 
b 2 — > from the analytical continuation of the measure to D < 2. The second is also an 
integrable singularity for b — > 0. They are handled by the following variable transforma- 
tions: 

/ f(l,b,c) = — — -%± f daa^sin^f- 2 

Jb Jy ' ' ; D- 7 D-12 S D Jo yi J 

[ dub~<f(i,b,c)e(i-b)e(c-b), (3.3i) 

Jo 
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D 

Figure 3.4: Numerical results for the diagram (3.29) 



D 


residue 


1.00 





1.05 


-0.6993 


1.10 


-1.5336 


1.15 


-2.6599 


1.20 


-4.4451 


1.25 


-8.0963 


1.275 


-12.111 


1.30 


-21.9321 


1.31 


-31.6537 


1.32 


-55.845 



where 



7T 



B = -aD-i 

1 2 



5D 



7 = u(d c - 2) = — - 2 
c = ^Jb 2 + 1 - 26cos(/3) 



(3.32) 

(3.33) 
(3.34) 

(3.35) 



The transformations are constructed in order to generate factors which compensate the 
singularities of the integrand. The factor b 1 in (3.31) e.g. exactly cancels the singularity 
6~ 7 in (3.29) as can be seen from the small b expansion (a = 1) of the integrand: 



4 + 2 
16£> 



(3.36) 



The parametrization with the angle f3 is convenient to disentangle the divergence for 
b 2 — > from the divergence for b — > 0. 

An additional difficulty arises as for small b (3.29) is the difference of two diverging 
terms. One can use the small b expansion, equation (3.36), in a domain determined by 
the program of approximately b < 10~ 7 . 

This integral is now performed by a simple numerical integration routine using Simp- 
son's rule. On a workstation, this integration takes some minutes for a precision of 10~ 4 . 

The result of the calculation can be found in figure 3.4. For D — > 0, the diagram 
vanishes, whereas for D — > 4/3 it diverges. This is due to the fact that for D > 4/3 a 
new subdivergence for & — > occurs. 
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3.6 The limit D 1 



The case D = 1 can be treated analytically. For D = 1 and omitting the counterterms, 
(3.29) reduces to 



+ 



4 Jo 
3T 2 (1 
~4 T(2 
for 4 



db(b-b 2 )- d < /2 
4/2) 



4 

: 3 



(3.37) 



The factor 3 in (3.37) is due to symmetry. The integral runs over one of three possible and 
equivalent sectors only. For simplicity the counterterms for the relevant subdivergences 
were omitted, as the integrand can by calculated analytically without counterterms and 
since they were constructed from a finite part prescription and thus give no contribution. 
This is indeed the the integral of the counterterms is: 



j f db (b- d -' 2 + (i - i,) 1 - 1 -' 2 + b- d -"(i - b) 1 - 1 '' 2 ) 



- ~ 

4 I 2 



+ 



4 ' 4 
for 4 = 3 



4 



+ 



r(f -4/2)r(2 -4/2)' 



r(3 - 4) 

These calculations show that the limit D — > 1 is correctly reproduced [21]. 



(3.38) 



4 Conformal Mapping of the Sectors 

In the previous section, we had to calculate some integral over distances, restricted to 
various sectors, such that one of the distances is larger than all the others and is set 
to L. In that case, the integrand was symmetric with respect to the three distances a, 
b and c and each of the three sectors gave the same contribution. In the following we 
shall deal with non-symmetric integrands and with integrals over more distances. We now 
introduce an extremely useful tool, the mapping of sectors, which will reveal its full power 
in the analysis of the 4-point divergences. Besides calculational convenience it also shows 
that the simple pole at 1-loop order is an universal quantity. We shall discuss this later. 
First we show how different sectors, arising in the analysis of the 3-point divergences, 
can be mapped onto each other. Let us remind that, with the measure (3.23) and the 
normalization (2.36), we had to compute an integral over some domain in the upper half 
plane, with the measure 

/ = %± r d Vl rdy 2 (y 2 f-\ (4.1) 

Jy J£) J — oo JO 

of a function / of the three distances a, b and c between the points F = (0, 0), E = (-L, 0) 
and G = (j/i, 2/2), given explicitly by (cf. figure 4.1) 

a = Z fixed, b = ^/( yi y + (y 2 y , c = J(L - j/x) 2 + (j/ 2 ) 2 . (4.2) 
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Let us consider that / is homogenous, with degree A, but not necessarily symmetric 

f(na } nb } kc) = k~ x f(a } b } c) . (4-3) 

The exponent A is called the conformal weight of the integrand /. 

The upper half-plane in y can be divided into three sectors A, B and C. The sector A 




Figure 4.1: The sectors A, B and C 

consists of all triangles A(a, 6, c) with a = L and 6, c < a as indicated in figure 4.1. This 
is the sector over which the integration was performed in section 3.3. The sectors B and 
C are the domains, where b and c respectively are the largest distances. 
We can map B onto A. This mapping consists of 2 steps: 

• the rescaling with respect to F by a factor L/\b\ which maps G onto G' and E onto 



E'. (4.3) implies that / is changed by a factor y-^ 

• a mirror operation, which maps G' onto E and E' onto E", leaving invariant F. This 
operation is a permutation of the first two arguments of /. 

The mapping G— s-E" is a special conformal transformation, the inversion with respect 
to the circle Sl(F). In complex coordinates it is 



V = 2/i + m 



y = yi + m 



y 



(4.4) 



One easily checks that this mapping B — > A is one to one. The measure in (4.1) transforms 
as 



d D y = d D y fy^j 



2D 



The final result i 



is: 



/ f(a = L,b,c)=f 



X-2D 



eA\L 



f(b,a = L, 



(4.5) 



(4.6) 



For the integrals which give the residue at the critical dimension, (4.2) is dimensionless 
and the degree A in (4.3) is A = 2D so that we have 



yes 



f(a = L,b,c) 



yeA 



f(b,a = L,c) . 



(4.7) 
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An analogous transformation is valid for the mapping of C onto A: 

f f(a = L,b,c)=f f(c,b,a = L) (4.8) 

We call k = X — 2D the conformal dimension of the integral. For integrals with conformal 
dimension zero, a conformal change of coordinates to map the various sectors simply 
permutes the vertices of the triangle. 

One word should be said about the conformal mapping for integrals which have not 
conformal dimension 0. As a is fixed to equal L one can always multiply the integrand by 
a power of a/L, by this way adjusting the conformal dimension to 0. Then the conformal 
mapping again consists in a pure permutation of the arguments. 

A more general method to look at the mapping of sectors consists in using the measure 
(3.7) of [18] over the distances, considered as independent variables (we set L = 1). The 
integral of / over (4.1) is 

f(a,b,c) = dfi D (a,b,c)x A {a,b,c)6(a-l)f(a,b,c) (4.9) 
yeA J 



with the measure dfiu defined as 



1 S D-1 j 2 7,2 -i 2 /oa/ L _^D-3 



dfi D (a,b,c) = ^da 2 db 2 dc 2 (2A(a, 6, c)f~ a . (4.10) 

8 JD 

A(a, 6, c) is the area of the triangle with edge lengths a, b and c and \ A the characteristic 
function of the sector A: 

X A (a, 6, c) = 9(a - b)9(a - c) (4.11) 

Substituting a = db } c = cb and using the homogeneity of / and of the measure, the 
r.h.s. of (4.9) becomes 

/ dd 2 db 2 drc 2 6(db - 1) b 2D ~ 2 ~ x (2A (5, 1, £)) D " 3 / (a, 1, c) x Aa, 1, c) (4.12) 
8 S " ' 



'D 



Substituting further b = b/d and using the fact that the ^-distribution restricts b to equal 
1 yields: 



/ dd 2 db 2 drc 2 6(1 - 1) (T 2D+A (2A(a, b, t)) D - 3 f(d } b, c) x Ad, b, c) (4.13) 
8 jd 1 



In order to get the same formula as (4.7), one considers functions / with conformal weight 
2D and renames b = a, a = b and c = c, to obtain: 

dfi£)(a, 6, c) 6(a — 1) /(a, 6, c) X^( a ? ^ c ) = y ^d(«, c) <^(« — 1) f(b, a, c) x B ( a ? ^? c ) 

(4.14) 

The result is equivalent to (4.8). 
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The method extends to the contraction of four or more points. We give the general 
result for N points connected by the distances x = Z\ } y = z 2} z 3} . . . , zn(n-i)/2 in a 
slightly different formulation: 

J d/i D (x } y } z 3} ... } z N{N _ 1) / 2 )f(x, y } z 3} ... } z N ( N _±y 2 ) S(x - 1) (4.15) 
= J d/i D (x } y } z 3} ... } z N{N _ 1) / 2 )f(x, y } z 3} ... } z N{N _ 1)/2 ) 6(y - l)x" (JV_1)D+A 

The ordering of the distances on both sides of (4.15) has to coincide. 

In the case of two or more contraction domains the same rules apply: one adjusts 
the conformal dimension of the integral to equal and then can freely choose the length 
which has to be fixed. The ordering of the distances of course has to be respected. 

The mapping of sectors immediately proves the universality of the pole term in e at 1- 
loop order, i.e. as long as no double pole appears, once the general normalization has been 
fixed: Demanding the longest, the shortest or any intermediate distance to equal L results 
in the same pole term. From these regularization prescription every other prescription can 
be constructed. This observation completes the proof. (The question whether different 
prescriptions to subtract the relevant divergences change the pole term is discussed in 
appendix A). 

5 3-body Interaction Renormalization: Easy Graphs 
5.1 Graph reducible to 2-point integrals 



a 



Figure 5.1: Distances in (5.1) 

The hrst contribution to the renormalization of the 3-point interaction comes from the 
MOPE (2.15). The divergent integral is: 



1 



a<L Jb<L Jc<L \a 2u c 2u + c 2u b 2u + b 2p 'a 21 



d/2 



(5.1; 



By the same procedure as in section 3.5, the pole term is extracted by fixing one by one 
the distances a, b and c to equal L and to be the largest. Using the mapping of section 
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Figure 5.2: Numerical results for the diagram (5.3) 

4 this integral is converted to an integral, where a = L = 1 fixed, b and c now running 
from to oo. The pole term is obtained by setting d = d c : 





1 



c 2u _|_ c 2vft2v _|_ yiu 



d c /2 



(5.2) 



The integration can easily be done as the integrals over both b and c simply yield Beta- 
functions: 



(2 - DY r 

v ' L \2 2-D 



(5.3) 



For D = 1 the integral is 2tt and decays rapidly for D — > 2, where it vanishes. A plot and 
some numerical values can be found in figure 5.2. 
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5.2 Graph reducible to a 3-point integral 




Figure 5.3: The distances in (5.4) 

The next diagram, which has to be calculated, comes from the MOPE (2.16). The diver- 
gent integral is: 




z<L Je<L 



d/2 



Ab 2u e 2u + (a u + c v + b v ) (a v + c v - b v ) (a v - c v + b v ) (c v + b v - a v ) / 



(5.4) 



The pole term is again derived by applying the operator L-^\d=d c to the diagram and 
leads to look at the sectors where a, 6, c or e equals L and all other distances are smaller. 
By mapping the sectors, it can be transformed to an integral over the sectors, where the 
largest of the distances a, b or c equals L, e now running from to oo. 

For the relevant divergence for 6 — > 0, the counterterms are found by applying a finite 
part prescription. This gives: 





+ / + 



--b- vdc 

2 



a + e + c + e 



The integral over e can still be performed analytically: 



i r 



r 



D 



D 

2-D J x V 2(2-1?) 



2-d r f 3D 

1 V2(2-D) 



-D 



+ + 

L=a>6,c J L=b>a,c «/L=c>a,i 



(a v + c v + b v ) {a v + c v - b v ) {a v - c v + b v ) {c v + b v - a/) / 

2^ 



d c /2 



\b 2v e 2v + (a v + c v + 6^) (a v + - 6^) (a v - c v + 6^) (c u + b v - a v ) 



D 



2(2-D) 



2 V 



(5.5) 



(5.6) 
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Figure 5.4: Numerical results for the diagram (5.7) 
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2.4618 


1.20 


2.4854 
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2.8518 


1.275 


3.4084 


1.30 


4.9365 


1.31 


6.5140 


1.32 


10.5032 


1.325 
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This integral is numerically integrated in its symmetrized version: 

i rf^lr^ D 



2 - D 



2-D 

r 



2(2- D) 



3D 



2(2-D) 



L=a>c>6 

4 



2(a- D + b- u + c 



D 



-D 



D 



2(2-D) 



\(a v + c v + b v )(a v + c v -b v )(a v -c v + b v )(c v + b v -a v ) J 

-a~ 3D / 2 (c~ D / 2 + b~ D / 2 ) - b~ 3D / 2 (a~ D l 2 + c- D > 2 ) - c~ 3D / 2 (b~ D l 2 + a- D ' 2 ) (5.7) 

The numerical problems and their solution are similar to those discussed in section 3.5. 
The singularity for a = 1 and 6 — > is: 



d r 



D 
2-D 



r 



D 



4(2 - D) 



r 



3D 



2j2zEH b 2-5D/2 ^ + 0^1-2^ + Q ^ 



(5.8) 



,2(2-D 

This equation determines the exponent 7 for the transformation of the measure (3.31 

5D 



7 



(5.9) 



as before. It also serves to eliminate the small b difficulties as discussed at the end of 
section 3.5. 

The explicit numerical results are given in figure 5.4. For D — > 1 the numerical value 
approaches 7r, which is the result obtained by analytical calculation. For D — > 4/3 it 
diverges, again reflecting the fact that for D > 4/3 a new non-integrable divergence 
appears for small b. 
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6 3-body Interaction Renormalization: Hard Graph 
6.1 The 4-point diagram 




Figure 6.1: The distances in (6.1) 



The last diagram, which contributes to the renormalization of the 3-point interaction, is 
depicted in figure 6.1. The associated divergent integral is: 



2v r2v 



1 



,2v 2v j2v , 12v 



b 2l/ + d 



-d/2 



. - .. (6 - 1} 

Applying the operator L-^j- from section 3.3 to extract the pole term gives 6 different 
sectors: Each of the distances a, 6, c, d, e or / might be the largest. In order to simplify 
the calculations, the mapping of section 4 shall be used. We want to end up with an 
integral, where c = L = 1 is fixed, whereas the other distances may vary freely: 




2is r2u 



1 



2u 2v t2u i i2u 



b lV + d l 



-d c /2 



(6.2) 



For c = 1 the measure, given by (3.25) and (2.36), simplifies to an integral over the vectors 
a and /: 



Sp-lSp-2 



da\ I da 2 a 



D-2 



dfi I df 2 I df 3 f°- 3 F(a,b,c,d,e,f) (6.3) 



We herein abbreviated the integrand by F(a, 6, c, d } e, /). 

For small a or /, F possesses relevant subdivergences. They are eliminated via a finite 
part prescription, if one uses the modified 3-body interaction (3.13): 



2is r2u 



1 



: 2u _ e 2u _ b 2u + d 2, 



-dc/2 



-vd c 



(6.4) 



The integrand now is integrable. 
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6.2 Improvement of the measure 



For D < 2 the measure defined in (6.3) is a distribution and suffers from a relevant 
divergence for f 3 — > 0. Geometrically these are configurations, where the tetrahedron 
spanned by a, 6, . . . , / has volume 0, i.e. is restricted to a plane. A finite part prescription 
has to be applied in order to make the measure finite. 

One may think of implementing this prescription by subtracting the singularity. This 
method however imposes at least numerical difficulties. It is better to eliminate the 
singularity by a partial integration with respect to / 3 , which is mathematically equivalent. 
As only d, e and / depend on / 3 , the integral 



df 3 f°- 3 F(a,b,c,d,e,f) (6.5) 



can be converted to 



= 7^ I df 3 /i^RFK b, c, d, e, f ) (6.6) 

where R is defined via 

r> 1 9 1 9 1 9 i n\ 

R -ddd + ~e¥e + Jdf (6 - ?) 

The strength of the divergences for a — > or / — > is unchanged. It is important to 
remark that this trick can not be used to eliminate the relevant divergences when a — > 
or / — > 0. It works for the divergence in / 3 , because the integrand does not directly 
depend on f 3 but on d } e and /, which themselves depend on f 3 . So the derivation of F 
with respect to f 3 does not produce a factor l// 3 but factors 1/d, 1/e and I//, which are 
not singular for f 3 — > 0. Explicitly: 



yd 

RF(a,b,c,d,e,f) = — 



1/2 (d~ D - e~ D ) (> - e 2v - b 2v + d 2v ) - a 2v f 2v 
a 2v j 2v - l - (c 2v - e 2v - b 2v + <P) 2 ) ^ 1 + a~ vd f- vd 



(6.8) 



6.3 Parametrization of the measure 

The main singularities for small distances appear for a or / small. As in section 3.5 
we therefore want to parametrize the measure with the help of these distances. The 
divergences for small volume of the tetrahedron spanned by a, . . . , /, a 2 — > or f 3 — > 
shall be treated by a parametrization in angles as by this way small distance and small 
volume singularities are best disentangled. We have chosen the parametrization indicated 
in the figure. One triangle is spanned by c and a with an angle f3 between, another by 
c and /, where the corresponding angle is a. The angle between the planes spanned by 
these two triangles is r. The distances as functions of a, / and /3, <r, r are: 

b = \j a 2 + I — 2a cos (3 
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c=l 

Figure 6.2: Parametrization of the tetrahedron 



e = Jp + 1 -2/ cos a (6.9) 



d = y (a cos /3 — 1 + / cos cr) 2 + (a sin f3 — f sin cr cos r) 2 + (/ sin r sin cr) 2 

The integrals over a and / run from to oo, the integrals over /3, a and r over the interval 
[0,7r]. As we do not want to map all the points which are far away, we have to find a 
reparametrization of the measure which behaves for a — > like a 1 and for a — > oo like a w , 
by this way eliminating the principle divergences. If u is equally distributed we can use 

i l 

a = u D -f(l - u)d-u . (6.10) 
The integral over f3 will be parametrized as 

{7T 

-(2a)D-i a < 0.5 
2 vr _L_ (6-H) 

7T - -(2 - 2a)D-i a > 0.5 

which differs from (3.31) and following by the different range of integration: [0, it] instead 
of [0,7r/2]. 

The integral over / (remind the factor /| came from the partial integration with 
respect to / 3 ) 

/ dhdpdPf^i fl) (6.12) 

can be written as 

" dff ' 1 I da{sma) D I dr (sin rf-^f) . (6.13) 



l 



Again it will be transformed 



/ = v D -~f(l - v)d-w (6.14) 
Furthermore we choose in the same spirit as for f3 

-(2rj)D+i 77 < 0.5 

Tr - -(2 -2rj)D+i n > 0.5 
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and 




C < 0.5 
C > 0.5 



(6.16) 



So the complete integral over four points is: 



Sd-iSd- 



7T 



SI (2-D)(D-l)D(D + l) 



1 1 1 
+ 



1 



du a D , 
o \ D — 7 u to — D I — u J Jo 



2-D 



1 1 

— + 



1 



1 



o \D-7t) u — Dl — vlJo 



da min(2a, 2 - 2a) 0-1 (sin /3) 

—D 

drj min(2?y, 2 — 2rj) D + x (sin cr 



D-2 



l-D 



JCmin(2C,2-2C) o (sin r)" D_1 / 2 R_F(a, 6, 1, J, e, /) 



(6.17) 



Another way of parametrizing consists in replacing the integral over the vectors a and 
/ by the integral over the vectors a and d. This parametrization is especially useful to 
eliminate divergences, when a and / simultaneously go to infinity. It will be used for the 
integration of one of the sectors in the next section. The new formulas are given here, a 
prime indicating new angles and distances as can be deduced from the figure, a' and r' 




c=l 

Figure 6.3: Alternative parametrization of the tetrahedron 

obey the same relations as a and r. As before /, now 

1 1 



d' = v D -i(l - v)D- 

The other new distances are: 

e 1 = \/ a 1 + d' 2 — 2ad' cos a 1 



(6.18) 



(6.19) 



f = y (d' sin a' sin r') 2 + (d' cos a' — a + cos /3) 2 + (sin f3 — d' sin a' cos r') 2 (6.20) 
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For the integrand (6.4), the exponents are: 

7 = 



to 



D 
W 



(6.21) 
(6.22) 



6.4 Decomposition into sectors 



Although the measure absorbs the principal singularities it cannot handle all of them. 
There remains e.g. a singularity for b — > and e — > 0. Two methods may be applied 
to handle the remaining integrable singularities. The first consists in using the second 
measure of section 6.3. The second is to map again some parts of the domain of integration. 
Thereby, we face the problem that the measure is no longer symmetric in the distances, 
as we have changed it in order to eliminate the relevant singularity for f 3 — > 0. We will 
therefore modify the measure (6.17) to 



Sd-iSd- 



7T 



s 2 



D 



(2 - D)(D- 1)D(D + 1) 



du a 



D+2 



1 1 



D 



7 u 



+ 



1 



1 



dvf 



D+2 



D 



1 

— + 



to — D 1 — u / 
1 1 N 



2-D 



7 v to 

l-D 



D 1 



v 



da min(2a, 2 - 2a)D-i ( sm /3) 

-D 

drj min(2?7, 2 — 2rj) D+ 1 (sin a 



D 



d( min(2C,2 - 2() d (sinr 



(6.23) 



which is, except for the geometric prefactor, the invariant measure in D + 2 dimensions. 
The integrand accordingly has to be modified to 



T 



1 



{2A(a,b,c)y 



KF 



(6.24) 



where A(a, 6, c) = |ac sin(/3) = \\J{a + b + c)(a + b — c)(b + c — a)(c + a — b) is the area 
of the triangle spanned by a, b and c. In our case 



T 



(a + b + c)(a + b- c)(b + c - a)(c + a - b)p 
1 



- 2 f (*-° 



-D 



,2i/ 



2v r2v 



1 

4 IC " 



,2v 



„2v 



b 2 " + d 



2v 



b 2v + J 2 ^) - a 2 ^ 2 ") 

dc/2-l 



-vi c 



(6.25) 



The integrand now is conformal invariant, as follows directly from equation (4.15). 
The sectors are decomposed as follows: 

1. (a < 2 or / < 2) and (b > \ or e > \) 

This sector is convergent: F(a, 6, c, d } e, /) is integrated directly, using the simple 
measure (6.17). 
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2. a > 2 and / > 2 

The measure (6.17) does not eliminate the singularity, when both a and / simul- 
taneously go to oo. The easiest way to integrate this sector is to use the second 
measure (6.18) ff. of section 6.3. 

The divergences of the integrand could also be eliminated by a mapping. This 
however induces new singularities due to the measure (the term l/(2A(a, 6, c)) 2 in 
T, equation (6.24)). This would not be the case, if we had not been forced to use 
the trick of integrating by parts the measure. 

3. b < | and e < | 

In this sector the mapping can be used successfully: a has to be exchanged with b 

and e with /. We get 

T(b, a, c, d, /, e), a < | and / < |, 

using the measure (6.23) 



6.5 Numerical integration 

6.5.1 General remarks 

In this section we want to discuss the algorithm used to perform the numerical integration. 
Let us recall that we have to treat integrals like (6.4) with the measure (6.17) and some 
modifications due to mappings of the sectors. These integrals run over the unit- volume 
[0, l] n with n = 5. Such high dimensional integrals cannot be done directly by means of 
Simpson's rule or derived methods as the integration time diverges like (1 / 1 precision) 71 " , 
a < 1 and a = 1 if no special analytical behavior is given. Normally Monte Carlo (MC) 
integration is a good choice, as the integration time grows like (1 / 'precision) 2 independent 
of the dimension n. However in its original form it is not applicable to our integral since 
in large domains of integration the integrand is nearly and in small domains much larger 
than 100, always assuming that the integral is normalized to 1. Therefore we implemented 
an adaptative Monte Carlo integration (AMC) procedure suggested to us by J. M. Drouffe 

[22]. \ 

The algorithm will be described schematically for n = 2: 



Figure 6.4: Schematical description of the AMC-algorithm 
We have an integrand f(x) where the main contribution to the integral 



I 



d n x /(; 



(6.26) 
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comes from a small domain, in which the integrand is large. In the figure this is indicated 
through a high density of points. We now try to integrate the square [0, l] n whose area is 
A = 1 with, say, 100 MC-sample points. The routine returns us the first moment /, Af 

being an estimate of the integral and the connected second moment f 2 — f 2 , A\J f 2 — f 2 
being an estimate of the deviation of Af from the integral /, henceforth called error 
estimate. A large error estimate indicates that the integral has not been found with a 
sufficient precision due to the fact that the integral is located in a small domain. Therefore 
we divide the domain of integration in 2 n subdomains (subboxes) as indicated in the figure 
and perform again a MC-integration with the same number of sample points in every 
subbox. The estimate of the integral now is 

( 6 -27) 

i 

where the sum runs over all boxes and fi and A{ are the mean value of / in and A{ the 
area of the z'-th subbox respectively. The new error estimate for the z'-th subbox is: 

Sli = A^J 2 -] 2 (6.28) 
This determines the total (statistical) error estimate: 

thtat = ^SI 2 (6.29) 

This time only in the upper left subbox the error estimate will be large and this box has 
to be subdivided again. We repeat the procedure recursively until the error estimate of 
every subbox is smaller than a given e. In our example this process stops after 3 recursive 
subdivisions. 

It should be emphasized that we did not divide the error estimate 61 in a subbox by 
\J number of samplepoints as can be done for Gaussian distributions. This is due to the 
fact that we are normally far from the domain where the central limit theorem is valid 
(Gaussian domain). 

6.5.2 Capabilities and performances 

The algorithm is able to integrate singular integrands. In low dimensions (n = 1,2) this 
works with appropriate accuracy in a very short time (some seconds). In higher dimensions 
(e.g. n = 5) however singular integrals may impose severe problems: The algorithm 
can simply overlook divergences. Of course, one could start with a finer sublattice, i.e. 
demand a minimal number of subdivisions in every direction, or one can use a smaller e. 
In practical calculations however, the execution time will explode. 

As a general rule the algorithm seems to work well for functions whose second moment 
/ f 2 (x) exists. For such well-behaved integrands the integration time scales slightly better 
than inverse proportional to the demanded precision. 

6.5.3 Implementation of the algorithm 

The algorithm works recursively, so the appropriate language is C. It was implemented 
on a SUN-workstation. A portation to FORTRAN is possible, but tedious because of its 
missing capabilities to use structures and recursions. 
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6.5.4 Numerics for the 4-point integral 
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Figure 6.5: Numerical results for the diagram (6.2). The first error in the table is the statistical error 
(6.29), the second the systematic error (6.30). 



We recall that 




was decomposed into 3 sectors. The main 



numerical problem, identical in each sector, is the concentration of the integral in small 
domains. This problem is solved by the AMC-algorithm discussed above. Additional 
numerical problems arise for D < 1.3 and become dominant for D < 1.1. The floating 
point routines of the workstation are implemented in double-precision (64 Bits, 52 Bits 
fraction 15 digits). In this range the routines more and more often fail to calculate the 
integrand, i.e. R_F (6.8) times the measure term, compare (6.17). These problems reflect 
the fact that for D — > 1 the (modified) measure becomes a distibution. They either return 
a wrong result or a so-called nonalgebraic number (NaN), ±oo or — oo. Only the latter 
events can be excluded from the integration. They are counted and serve as an error 
estimate via the formula 

#{NaN,±oo} 



81 



sys 



IP 



^{integration points} 



(6.30) 
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where p was determined from similar, but analytically known integrals to be 

/>«2. (6.31) 

These estimates have to be taken with some precautions as p depends on the actual 
integral and on the machine. The systematic errors in figure 6.5 were calculated by this 
method and seem to give a reasonable estimate for D — > 1, the limiting case for which 
the analytical result is 2tt (see next section). 

This line of arguments could be confirmed by working in quadruple precision (128 Bits, 
112 Bits fraction 33 digits). The systematic error for D = 1.05 was reduced by a factor 
3. No further analysis however was undertaken as the performance of the workstation 
drops by a factor 300. A reasonable calculation which before took one hour now needs 
two weeks. 



6.6 D^l 



In the limit D — > 1, 




can again be calculated analytically. 



This calculation is interesting as it reveals the connection to standard polymer theory and 



the fact that 




decomposes into three topologically different 



and non-equivalent diagrams. As in section 6.1 we keep c = 1 fixed. By a direct calculation 
it can be verified that the measure indeed reduces to an integral over a line. On this line 
two points, the endpoints of c, are already fixed. Then there are 12 different possibilities to 
distribute the last two points. They still can be separated into 3 topological inequivalent 
classes A, B and C, cf. figure 6.6. These are the three standard diagrams, arising in 
polymer theory. In each of these classes the line with c = 1 may be chosen to be the line 
connecting (12), (14), (23) or (34). Readers more familiar with Feynman diagrams arising 
in the framework of a scalar held theory may recover the three corresponding diagrams 
after a de Gennes transformation [23]. They contribute to the renormalization of the 
(^-interaction at d = 3 and are represented on the r.h.s.. Diagrams in one class can be 
mapped onto each other by the now well-known mapping of sectors. 

As an example, we calculate the diagram A, where c = 1 is chosen to be the line (23). 
Counterterms are neglected: 



1 

4 Ji 



da 



db (ab — V 



-3/2 



da , 

— a- 1 



2/1- 



T(l) 



-IT 



-1/2 



(6.32) 



The other three diagrams of this class were explicitly checked to give -| too. In the classes 
B and C all diagrams equal 0. So together, the contribution is 2tt confirming the numerical 
results. 
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Figure 6.6: The 3 topological inequivalent classes A, B and C 



This result was already obtained in a different parametrization in [24]. 

7 Results for the 3-body Hamiltonian 
7.1 Renormalization and critical exponents 

In this chapter we detail the renormalization of the model at 1-loop order, the derivation 
of the renormalization group equations and the determination of the critical exponents, 
along the line of section 2.3 and section 2.4. This is nothing but an application of the 
general procedure of [15, 16]. The purpose of renormalization is the following. We start 
from the bare Hamiltonian (2.1) 

K M = 2^ / \^r (x)f + g jf jf jf S d (r (x) - r (y))S d (r (x) - r (z)) . (7.1) 

Expectation values of "bare" observables (9o[?"o], i.e. functionals of the bare held r , are 
defined as 

(Oo[ro]) go = [v[r ]0 [r ]e- n °J r o] . ^ = / v[t . q] ~K [r ] _ (? _ 2) 

Expressed as perturbative series in the bare coupling constant g , they are UV hnite for 
e > 0, where we recall that 

e = W - 2vd , (7.3) 

if the relevant UV divergences due to the 2-point operator are subtracted according to the 
prescription of section 3.2. They suffer from UV divergences for e = 0, which appear at 
1-loop order as single poles in e. Renormalization consists in reexpressing the theory in 
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terms of renormalized quantities - the renormalized field r and the renormalized coupling 
constant g - defined as 

r(x) = Z~ 1/,2 ro(x) 



g = Z- d Z- 1 ^g 



(7.4) 



(with fi the renormalization momentum scale) so that the expectation values of physical 
observables, expressed in terms of the renormalized held r and of the renormalized cou- 
pling constant g, are UV finite as e — > 0. In renormalized quantities, the bare Hamiltonian 
is rewritten as the renormalized Hamiltonian (see equation (2.21)) 



n 



R\ 



KM 



Z 



D 



1 



Vr(x)) 2 +gZ g f fff 



6 d ( 



r x 



r(y))6 d (r(x) — r(z)) 



and the expectation values of renormalized quantities are 

fV[r]0[r]e 



■KM 



(0[r})" = (0[r]) 



90 



-K[r] 



(7.5) 



(7.6) 



JV[r]e 

At f-loop order, i.e. to hrst order in g, the counterterms are, in the spirit of the minimal 
subtraction scheme, chosen to have pure poles in e of the form (2.28), (2.29): 



Z 



i + g- + 

e 



Z n 



i + g- + 

e 



(7.7) 



Expanding (7.6) in terms of the renormalized coupling constant g we obtain 
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(7.8) 



where we made use of the standard abbreviations for the expectation value in the free 
theory 



(O[r])o 



and for the connected correlators 



JV[r]0[r]e 2-D L\(^ r f 
JV[r]e~^D fx |( Vr ) 2 



(7.9) 



(AB) conn = (AB) - (A)(B) 
(ABC) conn = (ABC) - (A)(BC) - (B)(CA) - (C)(AB) + 2(A)(3)(C) . (7.10) 

The poles in e of the bare expectation values are given by the MOPE and we have for 
instance 



O 
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+ > »+ 



+ finite terms 



+ finite terms (7.11) 
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Similarly, the third term on the r.h.s. of (7.8) contains a single pole which comes from the 
contraction of two 3-point operators into one 3-point operator: 





There are three different kind of contractions. Including the number of inequivalent ways 
of contracting the vertices, one has 



6 x 




+ 36 




+ 9 



(7.12) 



This term also contains poles which come from the contractions of 3-point operators 

into the I-point operator -|- : — ^ anc ^ ^/^\^ r /^\ > — > ' ^ rs ^ P°^ es 

will be subtracted by the same wave-function renormalization which subtracts the I-loop 
divergence (7. II). The second pole corresponds to a 2-loop wave-function divergence and 
need not to be considered here. Thus we have 





+ higher order wave-function divergences 
+ finite terms 



(7.13) 



From (7.8), (7.11) and (7.13), one sees that the single poles cancel and that the renor- 
malized theory is finite at 1-loop order if the counterterms Z and Z g are given by (7.7) 
with 



b = 3 




(7.14) 



The renormalized observables satisfy the renormalization group equations. They are easily 
obtained by calculating the variation of the renormalized quantities with respect to the 
renormalization scale //, which corresponds to the scale at which the theory is probed, 
keeping the bare Hamiltonian, which represents the microscopic theory, fixed. The flow 
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of the coupling constant is given by the /3-function: 



/%) 



d_ 



-eg 

1 + 9-§7, lnZ g + d c9§- g ^Z 



go ~ d a 
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(7.15) 



with d c = 3D/(2 — D). The scaling dimension of the field, v(g), is similarly obtained from 
the flow of the dimensionless renormalized field 



ln(f) 
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go,r 

D a 
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go 



+ 9^ + 0(g 2 ) 



(7.16) 



v(g) is related to the "fractal dimension" dp(g) of the membrane through v(g) = D / 'dp(g). 

The analytical and numerical calculations of the previous sections show that the co- 
efficients a and b (7.14) are both strictly positive for 1 < D < 4/3. This implies that for 
positive (and at least small) e, the /3-function has a non-trivial IR-attractive fixed point 
at 

9 * = ^ + °("» <" 7 » 

This fixed point governs the large distance (small fi) behavior of the membrane at the 
©-point. The scaling dimension of the membrane at the ©-point is given by 
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fa 



Kg* 



D 



+ e 



2(ad c + b) 



+ 0(e 2 



(7.18) 



Finally let us recall that, while the /3-function and the anomalous dimension v(g) depend 
on the normalization of the Hamiltonian and on the choice of the renormalization scheme, 
the result for the anomalous dimension v* is universal. In particular, it does not depend 
on the normalization introduced in section 2.5. 



7.2 Results for the anomalous dimension at 1-loop order 

One possibility to make the expansion is to consider the dimension D of the membrane 
as fixed and to vary the dimension of the external space d around the critical dimension 
d c = 3D/(2 — D). This means to expand v* in e = d c — d: 

= , MD) + Om , M D) = <™, 

The value of the coefficient A(D), which determines the scaling dimension of the mem- 
brane at order e, is plotted in figure 7.1 for 1 < D < 4/3. 

In this region, A(D) is positive and the scaling dimension v* is larger than the scal- 
ing dimension of the free "phantom" membrane v = (2 — D)/2, as expected since the 
interaction should swell the membrane. The correction vanishes for D — > 1. This is in 
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Figure 7.1: Anomalous correction to i>* at 1-loop order (3-point interaction) 

agreement with the known results [21] for polymers at the ©-point: corrections to the 
scaling dimension hrst appear at two loops, i.e. at order e 2 . Interestingly, although indi- 
vidual diagrams diverge at D = 4/3, the correction for the anomalous dimension seems 
to be regular as D — > 4/3. We shall prove in the next section that this is indeed the case 
and that the correction can be calculated analytically at D = 4/3. 

8 L> = 4/3 

As we remarked in the calculation of some of the diagrams, new singularities appear at 

D = 4/3. They come from the fact that the modified 2-body operator • h •, defined 

by (2.41), which is irrelevant for D < 4/3 and e = 0, becomes marginal at D = 4/3 and is 
relevant for D > 4/3. As argued in [15], it is this last interaction term which is relevant 
to describe the point for D > 4/3 and d close to the critical dimension, which is now 
given by d' c = 2(W - 2)/(2 - D). 

In this section we hrst discuss the limit D — > 4/3 of the 1-loop results obtained for 
D < 4/3 by the 3-point interaction only and show how the modified 2-point interaction 
emerges from the 3-point interaction as D — > 4/3. This allows to prove that this limit is 
regular as far as physical quantities such as the anomalous dimensions are concerned. 

In the second subsection we briefly recall the 1-loop results for the modified 2-point 
interaction, which are a priori valid for D > 4/3. The limit D — > 4/3 also turns out to 
be regular and the 1-loop correction for the anomalous dimension can be continued to 
D < 4/3. 

Finally in the third subsection we study the domain around D = 4/3 and d = 6, which 
is characterized by a crossover between the 3-point interaction and the modified 2-point 
interaction, described by the full Hamiltonian (2.20). We show that at 1-loop order this 
case can be treated analytically and that the crossover between the two interactions is 
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understood. 



8.1 The limit D —> 4/3 for the 3-body interaction 



Regarding the analytical and numerical results for the integrals involved in the calculations 
for the 3-body interaction (figure 3.4, 5.2, 5.4 and 6.5), one remarks that the integrals 



-f- } and 





are diverging for D — > 4/3 whereas the two 



others, 



-|- } and 



two are expected to dominate in the limit D 



+ 

4/3. 



stay finite. Therefore the hrst 



Let us start with the integral 



-|- > , given by (3.17) and (3.9). This integral 



has a global divergence as e — > 0. The additional subdivergence as D — > 4/3 is due to the 



subcontraction of two out of the three points: 



This subcontraction 



is given by the MOPE (2.13). If we denote by (y,z) the pair of subcontracted points, it 
takes the explicit form: 



C({y,z}) y , 



c ({y> z }) = \y - z \ 



»x + D({y,z}) y » — H m x + ... 

D({y,z})= l -\y-z\-^ 



(8.1) 



We remind the notation of points and distances represented in figure 3.1. The hrst term 
is the most singular one and is subtracted by the finite part prescription described in 
section 3.2. The second gives the singularity at D = 4/3. Remind, that the factorization 
of the integrations is valid only if b = \y — z\ is (much) smaller than a = \x — y\. So an 
effective IR cutoff L 1 has to be introduced. It will be specified later. Integrating in the 
domain b = \y — z\ < L' , we get 



b=\y-z\<L' 




+ finite terms (8.2) 



where we used the notation, similar to that of section 3.3, for the singular coefficient 
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b=\y-z\<L' 
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4(2 - vd) 



D({y,z}) 



b<L' 4 



b -v(d-2) 



L 



/2-vd 



(8.3) 



which diverges as vd — 2 — > 0. Since we are interested in the pole given by the global 
divergence for e = 0, we can replace vd by 3D/2 and the pole term 1/(2 — vd) by 

2 1 . . . T 

It is this term, associated to the subcontraction — > • h • , which 



3 4/3 — D ' 

gives the dominant contribution as D — > 4/3. 
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In this limit, another interesting thing happens: (8.3) depends only logarithmically on 
L' . This solves the delicate problem of how to choose this regulator when one integrates 
over the other distances. One can simply take L 1 to be the distance a between the two 

points in • h • . This is natural since this distance is the only one available. It is 

equivalent to integrate in the sector b < a. 

In the integration over the remaining distance, the dominant contribution comes from 
the contraction • h • — > -f- , given by a MOPE of the form 



H({x,y})* + + 



H({x,y}) = ^\x-y\-^ , /({*, y}) = ~^\x - y\ 



\2D-2-vd 



(8.4) 



The integration over a = \x — y\ has to be performed for a < L. Since there are 3 
inequivalent sectors we end up with 
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3 / D({y, z}) I({x, y}) + less singular terms 

'b=\y — z\<Ca = \x — y\<.L 
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(8.5) 



This result agrees with the numerical predictions. 
A similar analysis shows that 
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To perform the limit D — > 4/3 it is convenient to reexpress g in (7.5) in terms of 

1 
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2(4/3 - D 
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(8.6) 



(8.7) 



and analogously for bare quantities. The renormalization group functions become in this 
limit 



/3(g) = fi 



d_ 



21 

-eg + —g 2 + O(f) 



90 

2 - D 1 



(8.8) 



(8.9) 



These functions are used to analytically calculate the function A(D) plotted in figure 7.1 
in the limit D — > 4/3. 

In addition an interesting and striking property can be remarked: g is the appropriate 
variable to analytically continue to D > 4/3. For e > the fixed point g* stays positive 
whereas in terms of the original coupling g the associated g* becomes negative. This fact 
will be further clarified in section 8.3. 
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8.2 Comparison with the modified 2-body Hamiltonian 
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Figure 8.1: Anomalous corrections to i>* at 1-loop order (A(D) from the 3-point interaction and B(D) 
from the modified 2-point interaction) 

Since the modified 2-point interaction • h • becomes relevant for D > 4/3, it is 

interesting to compare the 1-loop renormalization group calculations for this interaction 
with the results for the 3-point interaction. These calculations are very similar to the 
original calculation for the 2-point interaction and can be performed analytically at 1- 
loop order [25]. Here we recall the principle of the calculation and the results. 
Let us start from the bare Hamiltonian 
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2 - D 

The dimension of the coupling constant b is 

e' = 3D - 2 



(Vr (x)) 2 + bo (-A r )S d (r (x) - r (yj) 



x Jy 



vd 



(8.10) 



f8.ll] 



and the model has UV divergences, i.e. poles in e', for e' = 0, that is for d = d' c = 
2(3D — 2)/(2 — D). The model is made finite by defining a renormalized held r and a 
renormalized coupling constant b: 



r x 



bo 



8.12) 
8.13) 



At 1-loop order, the counterterms are found to be 

z= b 



z h 




(8.14) 
(8.15) 
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The residues 



+ 



an 



d 



are analytic functions of D for < 



D < 2, listed in appendix B. The /3-function and the scaling dimension vib) for the held 
r are: 



1 d_ 
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(8.16) 



The new /3-function has a non-trivial IR hxed point b* > for e' > and the scaling 
dimension of the membrane at the ©-point, z/*, can be expanded in e' = d' c — d for hxed 
D: 

2 - D 



+ e'B(D) + 0(e 



(8.17) 



The 1-loop coefficient B(D) has no singularity at D = 4/3. Its value for 1 < D < 4/3 
is shown in figure 8.1, where it is compared to the 1-loop coefficient A(D) from the 3-point 
interaction. At D = 4/3, the expansion is around the same critical dimension d c = 6 and 
it is interesting to note that, although the expansion parameters e and e' are different, 
the two coefficients A and B are very close, but not identical at D = 4/3. 



8.3 Mixing of 2- and 3-body interaction at D = 4/3 



and 



interchange their role. Below 4/3 



At D = 4/3 the two operators 



is more relevant, above 4/3 



marginal: We have the interesting situation of a system with two coupling constants. The 
associated bare Hamiltonian, one wants to renormalize, is: 



'. At D = 4/3 and d = 6 both operators are 
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x Jy 



6 d (r {x) - ro{y))S d (ro{x) - r (z)) 



. , (8-18) 

fx J y J z 

In contrast to standard perturbation theory, where in hrst order of the coupling constant 
the divergences are single poles, the leading singularity of ^ jy-^j ^ * s a double pole, 
due to the sequence of divergent contractions 



This prevents us from performing the renormalization in the standard way. Let us look 

at the problem from another point of view. As the modified 2-point interaction • h • 

renormalizes the elastic energy |(Vr) 2 , perturbations in this operator have to be controlled 
by a small coupling b as is done in the Hamiltonian (8.18). On the other hand, the 3- 



point interaction 



renormalizes the modified 2-point interaction 



via 



the 
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contraction 



• H • ) • H • 



Therefore there has to be a small parameter controlling the ratio of 



and • h • . 



These demands are satisfied by redefining the 3-point coupling constant g as b g . The 
bare Hamiltonian becomes 



/ \(Vro(x)) 2 + b I j y (-A r )6 d (r (x) - r (y)) 
+ b go I I I 6 d (r (x)-r (y))6 d (r (x)-r (z)) (8.19) 

J X J V J z 



Orders in b and g are counted equivalently. The canonical dimensions of the coupling 
constants b and g are: 



eb = 3D — 2 — vd , e g 



>d 



(8.20) 



As the perturbative expansion is valid in the vicinity of the critical point d c = 6 and 
D c = 4/3, the two small parameters e& and e g are independent. The theory will be finite 
in terms of the renormalized Hamiltonian 
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(8.211 



involving the renormalized quantities 

r(x) = Z~ x I 2 Tq(x) 

b = z-^ 2 ' 1 Z^ 1 /i- eb b (8.22) 
g = Z-WZ^/i—go 

where the renormalization factors have up to hrst order in b and g the form 

Z = 1 + p h q — 

with constants p and q. We draw the attention to the important point that pole terms 
in are always proportional to b as should be evident from dimensional arguments. The 
same is true for e g and g. It is also important to note that this would not be the case for 
the parametrization (8.18), nor if one there replaces g by what one might be tempted 
to do. 

In order to explicitly perform the calculations, we have as in section 7.1 to expand 
(0)^ h in g and b: 



o) + 



1 - z 

2 - D 



Zf,bfi 



o 



• H • 
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Z b Z g bgii 
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+ Z 2 Z g b 2 g/j 2 ^ 
+ higher order terms 



(8.23) 



The following diagrams contribute in hrst order of g and b at D = 4/3 (for more details 
cf. appendix B): 
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(8.24) 
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(8.27) 



They determine the renormalization factors at 1-loop order: 

£b 



Z n 



Z h 



3 g 7 b 

1 +7 + o- 

4 e 3 2 e 6 

3 o 6 
1 - -— + - 



4e 



£b 



(8.28) 
(8.29) 
(8.30) 



As usually we dehne the renormalization group /3-functions of the two couplings b and g 
as their variation with respect to the renormalization scale // at hxed bare parameters: 



Pb(b,g) = fi 
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bo, go 
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(8.32) 



Plugging in the dehnitions of b and g, we get two coupled linear equations in /3& and /3 3 
which can be solved after some algebra. They lead to the /3-functions at 1-loop order: 
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Pg(b,g) 



-bg + 5b 2 
3 



4 

11, 

- £ g9 + Y bg+ 4 9 



8.33) 
8.34) 



48 




The scaling function of the held v(b } g) becomes: 

1 d 

v(b. q) = v In Z 

= \ + \h (8.35) 

The system of equations (8.33), (8.34) determines four hxed points in the (g, b) plane. The 
physical couplings must correspond to a repulsive interaction at short distance, hence to 
the domain (b > 0, g > 0). One of the hxed points is IR-attractive, one IR-repulsive 
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and the other two have one attractive and one repulsive direction. For special values 
of the parameters £& and e g , fixed points may coincide. Passing through these special 
values describes the transition from one fixed point to another, resulting in an eventual 
non-analyticity of the critical exponent vib^g). We first list the different critical points 
visualized in figure 8.2. 

Pi: The Gaussian fixed point b c = and g c = 0: it is stable for £& < and e g < 0. 

P2'. The fixed point b c = and g c = |e 3 describes also a trivial theory, although g c has a 
nontrivial value. Indeed, regarding the action (8.19), we see that both interactions 
are renormalized to 0. Also the critical exponent v(b } g) equals that of the free 
(Gaussian) theory. The stability condition is e g > and £& + e g < 0. 

P 3 : The fixed point b c = |e;, and g c = 0: for this non-trivial fixed point only the modified 

2- point interaction plays a role. It is stable for £& > and 1 le& > I0e 3 . 

P 4 : The fixed point b c = ^-(e& + e g ) and g c = ^( — 1 le& + 10e 3 ) is the most interesting. 
Both couplings flow to a finite non-zero value. This point is stable for e& + e g > 
and I left < 10e g . It corresponds to the earlier derived fixed point for the case of a 

3- point interaction only in the limit D — > 4/3 from below. We will explain that in 
more detail below. 

Let us discuss the graphics of figure 8.2: We can distinguish 8 different regions in the 
(d } D) plane around the critical point (d c = 6,P* C = 4/3), named A to H. The separating 
lines are: 

1. s g = separating D,E and A,H 

2. + e g = between E,F and A,B 

3. £6 = separating F,G and B,C 

4. lleb = 10e 3 between C,D and G,H 

The flowgraphs in figure 8.2 correspond to these regions A to H, starting with region H 
in the upper left corner. 

The situation encountered in the discussion of section 8.1 corresponds to an expansion 
in d — 6 for D = 4/3 fixed. This direction lies in the sector H and is near to the transition 
line lleb = 10e 3 , separating the domain of attraction of the fixed points b c = |e&, g c = 
and b c = ^-(e& + e g ), g c = ^( — + lOe^). The correct value for the anomalous 
scaling dimension of the field is thus given by the calculation with the modified 2-point 
interaction only. The fixed point described by the calculations for a 3-point interaction 
only before taking the limit D — > 4/3 is equivalent to the non-trivial but unstable fixed 
point b c = ^-(e& + e g ), g c = ^( — + lOe^). This can be seen by calculating the flow 
equation for bg } which has the same fixed point as the flow equation (/3-function) for g 
in section 8.1, if b here and g there are identified. In this parametrization also v(b } g) 
here and v(g) there are equivalent. The fact that the direction of expansion is near to 
the transition line, at which the 2 fixed points coincide, explains why the anomalous 
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Figure 8.3: The renormalization flow of the expansion for constant D. The fixed points are marked 
with points. An additional circle indicates an attractive direction. The graph was obtained through a 
numerical study of (8.33) and (8.34). The length of the arrows is scaled as the square root of the speed 
of the renormalization flow. 

corrections for the scaling of the membrane differ so little. The flow graph representing 
this situation is drawn in figure 8.3. 

Coming back to the general situation depicted in figure 8.2, the flows are such that: 

• In regions C, D and E, the Gaussian fixed point Pi or the pseudo-Gaussian fixed 
point P 2 are IR stable. The modified 2-point and 3-point interactions are irrelevant 
and the large distance properties of the manifold at the ©-point are those of a free 
Gaussian manifold 

• In regions A, B and H, the fixed point P 3 , described by the modified 2-point in- 
teraction only, is IR stable. The 3-point interaction is irrelevant and the modified 
2-point Hamiltonian (8.10), discussed in section 8.2, is sufficient to describe the 
large distance properties of the manifold at the ©-point through an e^-expansion. 

• Finally, in regions F and G the fixed point P 4 , which contains a mixture of 3-point 
and modified 2-point interactions, is IR stable. As discussed above, this fixed point 
corresponds to the limit D — > 4/3 for the 3-point Hamiltonian (2.1). Therefore the 
pure 3-point Hamiltonian is sufficient to describe the ©-point in an e g expansion. 

If we extrapolate these 1-loop results, we obtain the picture already summarized in 
figure 1.1 for the ©-point as a function of the external dimension of space d and of the 
internal dimension of the membrane D: The (o?, D) plane is separated into three regions: 
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• For D < 2 and d large enough, both the 3-point interaction and the modified 2-point 
interaction are irrelevant. The ©-point is described by the Gaussian model. 



• For d < d c = 3-D/ (2 — D) and D small enough, the 3-point interaction is more 
relevant than the modified 2-point interaction and governs the ©-point. 

• For d < d' c = 2(3D — 2)/(2 — D) and D large enough, the modified 2-point interaction 
is more relevant than the 3-point interaction and governs the ©-point. 

At 1-loop order, the separatrix between these two domains is given by line number 4. 
(lie?, = 10e 3 , with e& and e g given by (8.20)), i.e. by the line 

d = 108 D — 138 (8.36) 

Thus, if we trust this picture far from the critical point (d = 6, D = 4/3), we expect 
that for 2-dimensional membranes (D = 2), the modified 2-point interaction will always 
be the most relevant to describe the ©-point, even for d < 6. One also checks that 
the modified 2-point interaction is less relevant than the standard 3-point interaction to 
describe polymers (D = 1) in two dimensions (d = 2) at the ©-point. 

Finally, let us stress that our analysis of the relevance of the two interaction terms 
leads to results drastically different from naive power counting or approximate schemes. 
Naive power counting predicts a separating line given by 

i = < 8 - 37 » 

and that for D = 2 the 3-body interaction is always more relevant than the modified 
2-body interaction. Flory-type arguments give a separatrix 

d = 3D + 2 (8.38) 

while a Gaussian variational approximation leads to 

d = 6 . (8.39) 

Both approximations predict that for D = 2 the 3-body interaction is relevant for low 
dimensions d (d < 8 and d < 6 respectively). 



9 Conclusions 

In this paper we discussed the renormalization of self-avoiding tethered membranes at 
the tricritical point. The 3-body repulsive interaction has been considered first. From a 
technical point of vue the calculations at 1-loop order are difficult and led us to the devel- 
opement of new technical tools. These tools have successfully been applied to membranes 
with intrinsic dimension 1 < D < 4/3 at the ©-point and gave the critical exponent v* 
at 1-loop order. They equally should apply for the 2-loop calculations in the case of pure 
self- avoidance, equation (1.1). 
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The limit D — > 4/3, d ~ 6 is especially interesting, since in this situation a crossover 
between the 3-body repulsive interaction and a modified 2-body interaction takes place. 
For D — > 4/3 we were able to analytically calculate v* . The expansion around the point 
D = 4/3 and d = 6 revealed new and interesting phenomena. Thanks to a new double 
e-expansion we were able to completely describe the structure of the renormalization 
group flow and of the fixed points. The crossover is understood as the passing from 
one IR attractive fixed point to another. This approach also settles the question which 
interaction is expected to be relevant to describe the ©-point for the "physical" systems: 
While for polymers it is the 3-body interaction, it is the modified 2-body interaction which 
is relevant for 2-dimensional tethered surfaces. 

Although the models considered in this paper seem to be rather complicated and 
cumbersome they possess a great mathematical and physical richness, which we hope will 
further be explored in the future. 
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Appendix 

A About the Finite Part Prescription 

In order to eliminate the relevant UV divergences which appear in subcontractions 
— •, we had in equation (3.13) defined a subtracted 3-point interaction 



by adding a counterterm proportional to the 2-point interaction. Another possibility is 
not to add such a counterterm, but to subtract the relevant divergences by a finite part 
prescription on the level of diagrams. Let us demonstrate that for the diagram involved 
in the elastic term renormalization (section 3). We start from expression (3.9) 

2 d (c 2 + a 2 -b 2 )b 2v + (a 2 + b 2 -c 2 )c 2v + (b 2 + c 2 - a 2 )a 2v 

B{{x,y,z}) = 



d 1 

2 



\(a v + b v + c v ) < K a v + b v - c v )(b v + c v - a i/ )(c i/ + a v - ¥ 

. (A - 1} 

which shall be integrated over the sector A, defined in section 4 and characterized by 
L = a > b } c. Divergences occur only for b — > and c — > 0. The finite part prescription 
amounts to: 

f.p. / B({x,y,z}) = I B({x,y,z})- I ±-(b^ d + c^ d )a D ^ d (A.2) 

J A J A JA,B,C 4:1) 
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In the counterterm the integration is not restricted by 6, c < a = L. Mapping the sectors 
B and C onto A yields a new expression c'(e) for the residue in (3.18). It differs from the 
one obtained in section 3, given in (3.21), by 

c'(e) = c(e) ~ jfij (a- vd b D - vd + c- vd b D - vd )(b- e - I) 

+ (a-' jd c D -' jd + b- vd c p - vd ){c- e - 1) (A. 3) 

This difference is of order e as e — > and so does not change the residue c(0) at e = 0. 

The same statement holds for the other diagrams involved in the 1-loop renormal- 
ization of the 3-point interaction. Different subtraction prescriptions of the relevant di- 
vergences (associated to the 2-point interaction) thus do not change the pole terms at 
e = 0, which determine the renormalization group functions at 1-loop order, but they 
may change the finite parts at that order. This implies that the renormalization group 
functions at 2-loop order will depend on the choice of the subtraction prescriptions. 

Another point, which has to be discussed, is whether the limit d — > d c in (3.21) or with 
the modified counterterm discussed above can be performed. It is clear that problems 
may only arise for small b or c. In this limit, the integrals may be expanded in the form 

/ b~ x (A.4) 

J b 

According to the definition of the finite part prescription it is sufficient to check that no 
pole term in 1/(D — A) may occur for D — A = 0. The reader may easily verify that this 
is indeed the case. 



B Some Diagrams 

In this appendix the MOPE is given for all contractions not calculated in the main text. 
In addition the residue of some more diagrams can be found, which was used but not 
calculated up to now. 

The labeling of distances is either unambiguous or follows the notations in figures 3.1, 
5.1, 5.3 or 6.1. 



d 



1 



2 \a 2v + b 21 



d/2 + l 



+ 



d + 2 / 1 



8 \a 2v + b 2 " 



d/2 



+ ... (b.i; 



3D -2 r 



D 
2-D 



3D-2-vd 



{2-Df r 



2D 
2-D 



:b.2) 



d + 2 r 



• H • 



D 
2-D 



2D-vd 



4(2 - D) 2 r (^L 



;b.3) 
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d{d + 2) ( 1 \ d ' 2+2 

a 2v _|_ tfu 



d 2 + 6d-8 ( 1 \ d/2+1 
+ h— • — h — . + 



16 



a 2u _|_ tfv 



10D — D 2 — 8 r 



2-D 



3D-2-vd 



2(2 -D) 3 r 



2D 
2-D 



d + 2 



2D-2-vd 



+ + ... 



+ 



d 



3D-2-vd 
I xd/2 + 1 



2 — D 



2 \a 2v + b 21 



w-2 r 



+ ... 



2-D 



3D-2-i/d 



(2-^) 2 r ^ 



;b.4) 

;b.5) 

;b. 6 ) 
;b.7) 

;b.s) 
;b. 9 ) 



d/2 



_ 46 2i/ e 2i/ + (a v + + 6^) (a v + - 6^) (a v - c v + 6^) (c v + b v - a v ) 



;b.io) 



• + 



■AD-2-vd 




■AD-2-vd 



;b.ii) 



C Renormalizat ion of Relevant Operators: Away from 
the Tricritical Point 

The tricritical point (©-point) was obtained by demanding the renormalized coupling t 
for the 2-point interaction to vanish. It is possible to study the general situation with 
t ^ by regarding the renormalized Hamiltonian: 



KM 



z 



2 — D Jx 2 



Uvr(x)) 2 +tV D+e)/2 / / S d (r(x) - r(y)) 

2 Jx Jy 

+gZ gt i s III 6 d (r(x) - r(y))6 d (r(x) - r(z)) (C.l) 

•J x J v J z 
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The dimensions of the renormalized couplings t and g are adjusted to vanish, g was 
defined in (7.4), t analogously is 



(C.2) 



From general arguments it is known that t as the coupling of a relevant operator does not 
appear in the renormalization factors Z, Z t and Z g as long as we stay within the frame of 
the minimal subtraction scheme. So Z and Z g are unchanged. In addition to (7.8) there 
is to hrst order in t: 

{0[r])l t ~ (0[r])f = -tZ tV PW (O * 

+t gi i DI2+3el2 (o 

This determines 

/ 1 • <i ( * (* 



+ 0(g 2 



(C.3) 



+ 0(g 2 



(C.4) 



The scaling of the 2-point interaction is thus described by the renormalization group 
7f-function: 



d 



to, go 
D + e 

D + e 



t 

-/%) 
+ 9R 



— — — In Z + — In Z t 

2 Og Og 



+ 0(g 2 ) (C.5) 



The new diagram involved in (C.5) has already been calculated: 





(C.6) 



This is understood from the fact that the additional exterior leg on the r.h.s. of (C.6) 
does not appear in the calculation of the MOPE coefficient. 

By now it should be clear that the limit D — > 4/3 is regular if one uses the rescaled 
coupling g in (8.7). It is equally possible to calculate around D = 4/3 and d = 6 following 
the discussion in section 8.3. The reader willing to perform this exercise will find the 
necessary diagrams in appendix B. 
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